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We report first-principles calculations of inelastic Seebeck coefficients in an aluminum monatomic 
junction. We compare the elastic and inelastic Seebeck coefficients with and without local heating. 
In the low temperature regime, the signature of normal modes in the profiles of the inelastic Seebeck 
effects is salient. The inelastic Seebeck effects are enhanced by the normal modes, and further 
magnified by local heating. In the high temperature regime, the inelastic Seebeck effects are weakly 
suppressed due to the quasi-ballistic transport. 



The electron- vibration interaction plays an important 
role in molecular electronics. Electrons flowing in nano- 
j unctions are characterized by quasi-ballistic electron 
transport Only a small fraction of electrons expe- 
rience the inelastic scattering. Electron-vibration inter- 
actions cause discontinuities in the current- voltage (I-V) 
characteristics known as the inelastic current tunneling 
spectroscopy (IETS) [2[. The IETS can provide informa- 
tion on the underlying atomic structures of junctions 
It also gives important signals to the molecular junction 
characterization |4|]. 

Electrons that travel with energies larger than the en- 
ergy of normal modes can excite corresponding vibrations 
in the nano-structure anchoring the electrodes. This 
effect causes local heating in the nano-structure d-Q- 
Heating occurs when electrons exchange energy with the 
excitation and relaxation of the energy levels of the vibra- 
tion of the nano-structured object that anchors the elec- 
trodes. The heating power is typically within 15% of the 
electric power (IVb) supplied by a battery even at ambi- 
ent temperatures because of the quasi-ballistic transport. 
The heat generated in the central wire region is dissipated 
to the bulk electrodes via phonon-phonon interactions. 
The heat generation eventually equilibrates the heat dis- 
sipation, where the wire region reaches an effective lo- 
cal temperature T w higher than the electrode tempera- 
ture T e . Local temperature depends on several factors: 
the strength of coupling between electrons and the vi- 
brations, the background temperature, and the thermal 
current which dissipates heat. 

In the last decade, remarkable progress has been 
achieved in measuring the Seebeck coefficients in nano- 
j unctions js _ U|- These experiments have shed light on 
the design of possible energy-conversion nano-devices, 
such as nanoscale refrigerators and power generators [12 1 - 
These experiments have also inspired rapid development 
in the theory of thermoelectric nano junctions [13l-l22|. In 
bulk systems, diffused electrons scattered by phonons can 
significantly affect the Seebeck coefficient. However, the 
effects of the quasi-ballistic electrons scattered by vibra- 



tions of the nano-structure on the Seebeck coefficient are 
relatively unexplored in nano junctions. To the best of 
our knowledge only one report for inelastic Seebeck coeffi- 
cients in molecular junctions based on model calculations 
is available in the literature [23|. In this Letter, we inves- 
tigate inelastic Seebeck coefficients from first-principles 
calculations. 

Following the work of Chen, Zwolak, and Di Ventra [Hj], 
the many-body Hamiltonian of the system is H = H e i + 
H V ih + H e i- V ib, where H e i is the electronic part of the 
Hamiltonian under adiabatic approximations and H V ih is 
the ionic part of the Hamiltonian, which can be casted 
into a set of independent simple harmonic oscillators via 
normal coordinates. The normal mode frequencies are 
(jjj, and R e i- vih is a part of the Hamiltonian for electron- 
vibration interactions which has the form of, 
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where a, /? = {L,R}; Mi is the mass of the z-th atom; 
Ai^j is a canonical transformation between normal and 
Cartesian coordinates satisfying Ai^jAi^j' = Sjj'', 

bj is the annihilation operator corresponding to the j- 
th normal mode, and a L ^ is the annihilation operator 
for electrons; the coupling constant J^ a ^ 2 between elec- 
trons and the vibration of the i-th atom in \i (= x, y, z) 
component can be calculated as, 

J&% =J dr l rfK||[^ lK|| (rr[^y^(r,R,)< K|| (r)] 
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where V ps (r, R^) is the pseudopotential representing 
the interaction between electrons and the i-th ion; 
^ < kk^ ,R \ t ) stands for the effective single-particle wave 
function of the entire system corresponding to inci- 
dent electrons propagated from the left (right) electrode. 
These wave functions are calculated iteratively until con- 
vergence and self-consistency are achieved in the frame- 
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FIG. 1: (color online) Feynman diagrams of the first- 
order electron-vibration scattering processes considered in 
this study. 



work of DFT combined with the Lippmann-Schwinger 
equation |24j . 



where a = {L,R}; S = (1) represents the process of 
phonon emission (absorption); the other two parameters 
in Eq. (jU) can be obtained by the relations Bfip) = 

—Bf'yp) an d ^^1(2) = ~-^ ] ji{2y the avera g e number of 
local phonons is (rij) = l/{exp[hujj/(kBT w )] — 1} , where 
T w is the effective wire temperature. 

The rate of energy absorbed (emitted) by the an- 
chored nano-structures due to incident electrons from the 
/3 = {T, R} electrode and scattered to the a = {T, R} 
electrode via a vibrational mode j is denoted by W^' 2 ^ . 
The total thermal power generated in the junction P, 
calculated from the Fermi golden rule, can be written as 
the sum of all the vibrational modes of eight scattering 
processes shown in Fig. [TJ 
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The rate of heat dissipated to electrodes via phonon- 
^honon interactions is calculated using the weak link 
model, 



where G is the Green's function of the biased bimetallic 
electrodes with Vb = (hr — where Hr^l) 18 the 

chemical potentials deep in the right (left) electrode, re- 
spectively; the wave function of the bimetallic junction, 
ek ( r )> * s calculated by solving a combination of the 
Poisson and Schrodinger equations until self-consistency 
is achieved, where the boundary conditions are given by 
the electrons deep inside the biased electrodes. The in- 
clusion of a single molecule bridging the bimetallic elec- 
trodes is considered as the scattering center, described 
by the potential V. 

Our starting point is the inelastic current when con- 
sidering electron- vibration interactions, 
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where K = 1.59 eV/ag is the stiffness of the 4-A1 atom 
chain connected to the electrodes obtained from the 
total energy calculation [25|; Nl(r){E) is the spectral 
density of local phonon DOS at the left (right) elec- 
trode surface from first-principles calculations [26|; and 
™l(r) = l/(e E/KBTL ^ - 1) is the Bose-Einstein distri- 
bution function. The effective local temperature T w is 
obtained when heat generation in the nano-structure and 
heat dissipation into the bulk electrodes reach balance. 

We calculate the inelastic Seebeck coefficient based on 
the inelastic current described in Eq. ((3]), which is a func- 
f E )~(B -B )]r(£), tion of Tl , Tr ^ and y B = ^ R _ We con sider 
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l/{exp[(E-ti L ( R) )/(k B T LiR) )] + l}b the 
Fermi-Dirac distribution function describing the statistic 
of electrons deep in the left (right) electrode with temper- 
ature Tl(r-) and chemical potential Hl(r)\ the transmis- 
sion function r{E) = J dR J dK^f^ - 
^^e*K\\ ^EK\\ ) * s calculated from the electronic part of 
the wave functions ^EK\\ • The terms B L ^ represent 
the corrections to the elastic current considering the eight 
first-order scattering processes depicted in Fig. [TJ 



an extra current induced by an infinitesimal tempera- 
ture difference (AT) across the junction. This current is 
counterbalanced by an extra current driven by a voltage 
(AV), which is induced by AT via the Seebeck effect, 
i.e., 
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After expanding the above equation to the first order in 
AT and AV, we obtain the inelastic Seebeck coefficient 
(defined as S e i+ V ib = AV/AT), 



where a, f3 = {T, R} and a ^ f3. The parameters B^ff, 
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and ^ denoted in Eq. (j4j) are, 
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FIG. 2: (color online) (a) Local temperature T w as a function 
of Vb for T e = 0, 4, 10, 50 K; (b) the differential conductance 
and the absolute value of dG/dV due the electron- vibration 
interaction without [solid (black) line] and with [dash (red) 
line] local heating as a function of bias for T e = 12 K. The 
schematic shows the normal modes that contribute to the 
jumps in the local temperature and inelastic current profiles. 
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FIG. 3: (color online) Elastic Seebeck coefficient [solid (black) 
line], inelastic Seebeck coefficient without local heating [dash 
(red) line], and that with local heating [triangle (blue) line] 
(a) as a function of bias Vb for T e = 4, 8, 12 K (in upper 
panel) and T e = 50 K (in lower panel); and (b) as a function 
of T e for Vb = K (in upper panel) and and Vb = 30 K (in 
lower panel.) 
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where a = {L,R} and ^f( 2 ) are given by Eq. (j5j). 
The other two terms in Eq. (|9]) can be calculated with 
the following relations ^f- = ^f-(L ^ i?) and = 

-g§-{L ^ i?), where L ^ R represents the interchange 
between R and L. We see that, in the absence of electron- 
phonon scattering, Eq. (|9)) recovers the elastic Seebeck 
coefficient described in Ref. (28| . 

We now apply our theory to investigating the inelastic 
Seebeck effects of four Al atoms bridging the bimetallic 
Al electrodes modeled as electron jellium with r s ~ 2. 
The 4-A1 junction is structurally and electronically sim- 
ple such that the first-principle calculations reported here 
can be performed with a high level of accuracy. It, there- 
fore, serves as an ideal testbed for comparing the predic- 
tions of theory and measurements in experiments. We 
compare the elastic and inelastic Seebeck coefficients as- 
suming that the left and right electrodes share the same 



temperature T e . In order to qualitatively show to what 
extent local heating affects the inelastic Seebeck coeffi- 
cient, we choose to display inelastic Seebeck coefficients 
with and without local heating. 

In the case of "without local heating", we mean that 
the heat generated in the wire region is perfectly dissi- 
pated to electrodes such that T w = T e . When including 
"local heating", the effective local wire temperature T w 
is higher than the electrode temperature T e . Fig. [2ja) 
shows T w as a function of the applied bias Vb for various 
T e . We note that three jumps occur at Vb = 2.5, 20, 
and 40 mV, corresponding to the energies of the normal 
modes. The sharp increase in T w at Vb = 20 mV corre- 
sponds to the first longitudinal vibrational mode. Two 
degenerate transverse modes are present in the x- and 
^/-directions at Vb =2.5 mV, and we show the represen- 
tative one in Fig. O Due to the selection rule, the con- 
tributions to local heating from modes with vibrational 
components perpendicular to the direction of electron 
transport (z-direction) are unimportant. For T e = 0, 
4, and 10 K, T w displays larger jumps at Vb = 20 mV, 
where eVs is the energy of the first longitudinal vibra- 
tional mode. For T e = 50 K, the signatures of normal 
modes in T w are wiped out by high temperatures. The 
increase in local temperature is less significant at higher 
T e . This is due to increasingly efficient heat dissipation 
caused by the increase of phonon population in the elec- 
trodes, as shown in Eq. ([7]). Fig.[2fb) shows the inelastic 
profile of the conductance (G = dl/dV) and derivative 
of conductance (d 2 I/dV 2 ) as a function of bias with and 
without local heating. Local heating enhances the ef- 
fects of the electron- vibration interactions on the inelas- 
tic current because of increased average number of local 
phonons. 



Figure [3^a) shows Seebeck coefficients as a function of 
the applied bias Vb for various T e . For each temper- 
ature, we calculate Seebeck coefficients in three cases: 
elastic Seebeck coefficients So, inelastic Seebeck coeffi- 
cients without local heating Si , and inelastic Seebeck co- 
efficients with local heating S2. The difference between 
the elastic and inelastic Seebeck effects is more salient 
in the low temperature regime around Vb = 20 mV 
[see the cases of 4, 8, and 12 K in the upper panel of 
Fig. 03(a)]. The profile of inelastic Seebeck coefficients 
vs. Vb displays a strong signature corresponding to the 
longitudinal vibrational mode at Vb = 20 mV, where 
the magnitude of the Seebeck coefficients are increased. 
This feature is related to the suppression of the inelas- 
tic current around Vb = 20 mV [Fig. Efb)], where the 
transmission function effectively decreases. This leads to 
larger magnitudes of Seebeck coefficients because S oc 

-t'( M )/t(m) & 

The inclusion of local heating enhances the effect of 
electron-vibration on Seebeck coefficients further. In the 
low temperature regime, the upper panel of Fig. [31(a) 
shows that Si (without local heating) significantly differs 
from S2 (with local heating). This is because of the large 
difference between T w and T e , as shown in Fig.[2fa). For 
Vb < 30 mV, T w and T e become almost identical when 
the T e is large. Consequently, the difference between Si 
and S2 becomes small [see cases of T e = 50 K in the 
lower panel of Fig. [31(a).] In all cases, the transverse 
modes are negligible to the inelastic Seebeck coefficients. 
Fig. 09(b) shows Seebeck coefficients as a function of T e for 
Vb =0 and 30 mV in three cases: So, Si, and £2. In the 
high temperature regime (T e > 50 K), the magnitudes 
of inelastic Seebeck coefficients (Si and S2) are slightly 
decreased compared with the elastic Seebeck coefficients 
(So) due to small probability of electron- vibration scat- 
tering. 

In summary, we investigated the elastic and inelastic 
Seebeck coefficients with and without local heating in 
the 4-A1 atomic junction using first-principles calcula- 
tions. In the low temperature regime, the signature of 
normal modes in the profiles of inelastic Seebeck effects 
is salient. The inelastic Seebeck effects are enhanced by 
electron- vibration interactions due to the drastic suppres- 
sion of the inelastic current at the bias corresponding to 
the normal mode with longitudinal vibrational character. 
Local heating enhances the inelastic Seebeck effects fur- 
ther due to increased average number of local phonons. 
In the high temperature regime, the inelastic Seebeck ef- 
fects are slightly suppressed by electron- vibration inter- 
actions due to quasi-ballistic electron transport in nano- 
junctions. The signature of normal modes in inelastic 
Seebeck coefficients and local temperatures is wiped out 
by the tail of the Fermi-Dirac distribution. 
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